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Abstract. We report on the detailed and systematic study of field-line 
twist and length distributions within magnetic flux ropes embedded in In¬ 
terplanetary Coronal Mass Ejections (ICMEs). The Grad-Shafranov recon¬ 
struction method is utilized together with a constant-twist nonlinear force- 
free (Gold-Hoyle) flux rope model to reveal the close relation between the 
field-line twist and length in cylindrical flux ropes, based on in-situ Wind 
spacecraft measurements. We show that the held-line twist distributions within 
interplanetary flux ropes are inconsistent with the Lundquist model. In par¬ 
ticular we utilize the unique measurements of magnetic field-line lengths within 
selected ICME events as provided by Kahler et al. [2011a] based on energetic 
electron burst observations at 1 AU and the associated type III radio emis¬ 
sions detected by the Wind spacecraft. These direct measurements are com¬ 
pared with our model calculations to help assess the flux-rope interpretation 
of the embedded magnetic structures. By using the different flux-rope mod¬ 
els, we show that the in-situ direct measurements of held-line lengths are con¬ 
sistent with a hux-rope structure with spiral held lines of constant and low 
twist, largely diherent from that of the Lundquist model, especially for rel¬ 
atively large-scale hux ropes. 
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1. Introduction 

Magnetic flux ropes are a type of well organized magnetic field structures embedded 
in space plasmas. The existence of such structures is best conhrmed by in-situ space¬ 
craft observations and the associated modeling when the structure is traversed by the 
spacecraft [e.g., Burlaga, 1995; Lepping et ai, 1990, 1997]. In addition, many studies on 
the origination of such structures also provide mostly indirect evidence to support such 
interpretation of these structures as magnetic flux ropes [e.g., Webb et ai, 2000; Long- 
cope et ai, 2007; Qiu et al, 2007; Demoulin, 2008; Qiu, 2009; Vourlidas, 2014]. They are 
found in Interplanetary Coronal Mass Ejections (ICMEs), the interplanetary counterparts 
of CMEs originating from the Sun. 

Some ICMEs are traditionally categorized as Magnetic Clouds (MCs) that possess a 
specihc set of signatures based on in-situ spacecraft measurements of both magnetic held 
and bulk plasma properties. A more modern view of all ICMEs containing hux ropes is 
also emerging [Gopalswamy et ai, 2013a, b; Xie et ai, 2013; Hu et ai, 2014]. This seems 
reasonable especially considering that most origination mechanisms for CMEs involve 
magnetic hux ropes no matter whether they are considered to be pre-existing prior to 
eruption, or generated during the process. Moreover the subsequent argument is that 
such structures originating from the Sun and propagating into the interplanetary space 
may not be properly detected by the in-situ spacecraft. Each spacecraft only provides a 
very localized, single-point measurements of the structure traversed. Therefore depending 
on the relative spacecraft path across the structure, the variability and limitation in the 
in-situ signatures of magnetic hux ropes are signihcant, resulting in the incidences when 
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the flux-rope structure is present, but the in-situ signatures are lacking [e.g., Jian et ai, 
2006]. However if one adheres to the traditional definitions of MCs, which satisfies these 
criteria: 1) relatively strong magnetic held magnitude, 2) smooth rotation in magnetic 
held direction, and 3) relatively low proton (3, the ratio between the plasma pressure 
and the magnetic pressure, one can likely derive a magnetic flux-rope structure from the 
in-situ data. 

Effort has been put on in-situ modeling of magnetic flux-rope structures in order to 
extend the current capability thus to better reveal and characterize these structures in 
a quantitative manner. Various flux-rope models utilize in-situ spacecraft measurements 
of magnetic field and plasma parameters along the spacecraft path and are based on 
either a cylindrical or toroidal geometry and magnetohydrostatic theory. They range 
from the well-known one-dimensional (ID) linear force-free cylindrical model [Lundquist, 
1950], to the corresponding toroidal model [Marubashi and Lepping, 2007; Romashets 
and Vandas, 2003], and to the fully two and a half dimensional (2|D) Grad-Shafranov 
(GS) reconstruction model [Hu and Sonnerup, 2002]. One particular model that has 
not been widely recognized is the so-called Gold-Hoyle (GH) model that was originally 
developed by Gold and Hoyle [1960] and was only applied in a limited number of studies 
[Farrugia et ai, 1999; Dasso et ai, 2006; Hu et ai, 2014]. The distinct features of this 
model, remaining ID, are that the field-line twist is constant across the radius and the 
corresponding equilibrium state is non-linear force free. 

In our latest study of Hu et al. [2014], we showed that the flux-rope structures as derived 
from the generally non-force free GS method are more consistent with the GH model than 
with the Lundquist model, especially in that the field-line twist distributions within IGME 
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flux ropes remain fairly constant for large-size, low-twist flux ropes. In the present study, 
we intend to elaborate more on this hnding and present additional consistency check by 
utilizing the unique measurements of held-hne lengths inside MCs. 

A unique set of in-situ spacecraft observations besides the magnetic field and plasma 
parameters in interplanetary space is the energetic electron burst onset. They appear 
as sudden increase in electron flux of energies up to a few hundred keV [Krucker et ai, 
1999; Kahler and Ragot, 2006; Wang et ai, 2011] as the electron beams propagate away 
from the source on the Sun to the location of the spacecraft along individual held lines 
connecting both ends. Under certain assumptions such as scatter-free propagation and 
coincidental release at the time of associated Type III radio burst, the path lengths of 
magnetic held lines can be derived especially inside MCs. There are two ways to obtain 
the length estimate based on electron burst onset observations: one is to directly calculate 
the length traveled by the product of the speed of electrons (of known energy) and the 
travel time (taken as the diherence between the onset time at 1 AU and the release time 
as given by the corresponding Type III onset time); the other is to linearly ht the onset 
times of electrons of diherent energies versus their inverse speeds (so-called inverse-beta 
method; [Kahler and Ragot, 2006]) and the slope yields the path length. The hrst study of 
comparing held-hne path lengths inside an MC utilizing the electron burst measurements 
was carried out by Larson et al. [1997]. They combined multiple in-situ measurements 
from the Wind spacecraft during an MC interval to derive held-hne lengths as measured 
by the energetic electrons travel time multiplied by the speed which were then compared 
with the lengths estimated based on certain hux-rope models of MCs. They found for one 
particular event that the path lengths at several locations inside the MC ranging from 
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about 1.2 AU near the center to about 3 AU near the boundary, consistent with flux-rope 
model estimates. Such type of study, rare but important, provided unique and direct 
supporting evidence for the interpretation of MC structures as magnetic flux ropes. Not 
until recently did Kahler et al. [2011a, b] extended that unique study by applying the same 
analysis to a set of Wind MC events and additional electron burst events from the ACE 
spacecraft. They derived the field-line lengths based on in-situ electron burst onset and 
associated Type III radio burst following the approach of Larson et al. [1997] and compared 
with two flux-rope models: one being the Lundquist model and the other flux-conservation 
model given in Larson et al. [1997]. Their comparison indicated poor correlation between 
the measured and the model field-line lengths with the latter being exceedingly larger, > 4 
AU with maxima reaching about 10 AU, especially for the Lundquist model. Their results 
cast doubt on the model fit to MC flux ropes by the Lundquist model which intrinsically 
possesses the property of increasing field-line twist thus length from the center towards the 
boundary of the flux rope at a rapid rate, approaching infinity at the boundary defined as a 
circular cylindrical surface of vanishing axial magnetic held. In addition, our own analysis 
[Hu et al., 2014] also showed that the field-line twist estimates from the GS method are 
not consistent with the Lundquist model but more aligned with the GH model of constant 
twist. In the present study, we will focus our analysis on the field-line length estimates 
based on the GS reconstruction results, supplemented by the corresponding estimates 
based on the GH model as well. 

Estimates of magnetic field-line lengths, by taking advantage of the unique and com¬ 
prehensive in-situ spacecraft measurements, not only provide constraint and validation of 
flux-rope models, but also provide possible measurement of one key parameter, the axial 
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length of a cylindrical flux rope. This parameter determines the quantitative measure¬ 
ments of the poloidal magnetic flux and the relative magnetic helicity contents [Qiu et ai, 
2007; Webb et ai, 2010]. Since all existing flux-rope models based on in-situ measurements 
are 2D at best in geometry, signihcant uncertainty exists for the axial dimension. An ef¬ 
fective axial length, L^ff, has to be used in order to determine the quantities of poloidal 
magnetic flux and relative magnetic helicity within a cylinder of hnite length Le// that are 
equivalent to the corresponding content contained within the actual flux-rope structure. 
In our effort to connect the ICME flux ropes with their solar source region properties, 
particularly by comparing the magnetic flux contents at both ends, a somewhat arbitrary 
range of Le// G [0.5, 2.0] AU was used [Qiu et ai, 2007; Hu et ai, 2014]. We strive to 
gain more insight and to obtain a rehned range of effective axial length from the current 
analysis of magnetic held-line length estimates inside MCs. 

The article is organized as follows. We present the detailed description of magnetic 
held-line length estimates inside the MCs in the next section, for the Wind spacecraft 
MC events given by Kahler et al. [2011a] for which the length measurements based on 
electron burst onset were published. We will reconstruct the structures of these events 
and derive the relevant characteristic parameters by using the GS method and the GH 
model will be primarily utilized to provide extrapolated estimates on held-line lengths. 
The approach of obtaining various length estimates is described in Section 2. These 
estimates are compared one-by-one with the corresponding electron burst measurements 
from the Wind spacecraft. Three cases are selected to be presented in detail in Section 3. 
A summary of our results and comparison for all events is given in Section 4. We hnally 
conclude and discuss the implications of our results in the last section. 


DRAFT 


February 19, 2015, l:21ain 


DRAFT 



X - 8 


HU ET AL.: MAGNETIC FIELD-LINE LENGTHS IN FLUX ROPES 


2. Magnetic Field-line Length Estimates 

We re-examine the events listed in Table 1 of Kahler et al. [2011a], total of 8 Wind 
spacecraft MC events with given electron burst measurements. We are able to successfully 
reconstruct 7 MC events by the GS method, except for the one on 2 May 1998. Therefore 
this event is excluded from our analysis. In addition, only the measurements of electron 
burst events occurring in the identified GS intervals (see Table 1) are utilized in our 
analysis. Others falling out of the GS intervals are excluded as well. 

2.1. GS Reconstruction Results 

The GS reconstruction is to solve the plane Grad-Shafranov (GS) equation of the mag¬ 
netic flux function A{x, y) (equivalent to the component of the magnetic vector po¬ 
tential; in direct analogy to a stream function in two-dimensional incompressible flow) 
on the cross-section of a cylindrical flux rope. Therefore the transverse magnetic held 
{Bx, By) is fully determined by the hux function A and the magnetic held lines are lying 
on cylindrical iso-surfaces of A, called A shells. The non-vanishing axial magnetic held 
component Bz becomes a function of A only, yielding a cylindrical but non-axisymmetric 
hux-rope conhguration. Detailed descriptions of the method and the latest updates were 
reported in Hu and Sonnerup [2002]; Hu et al. [2013, 2014]. 

Various physical quantities characterizing such a hux-rope structure can be derived 
including the axial magnetic held Bz, the axial electric current density and current, the 
toroidal (axial) and poloidal magnetic hux ^t,p, the relative magnetic helicity K^, and 
the magnetic held-line twist. They are all functions of A alone and vary across distinct 
A shells. Figure 1 shows the summary plots of these quantities as they vary along the A 
shells for all the Wind MG events we examined. They generally exhibit a similar pattern 
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to the other flux-rope events as we first reported in Hu et al. [2014] in this congregated 
manner. The magnetic fluxes increase monotonically from the center toward the outer 
boundary while the poloidal flux is generally larger than the toroidal flux. The relative 
magnetic helicity also increases monotonically and smoothly. So does the electric current 
since they are all accumulative integral quantities. The axial held, on the other hand, 
shows a monotonically declining profile from the center outwards, typical of a flux-rope 
structure. The maximum value (i^zo) ranges between a few and a few tens nT. The most 
irregular variation exists in the current density which represents the first-order derivative 
of a transverse pressure with respect to A. The field-line twist estimates displayed here, 
i.e., th = \Kr\/^f and tf = are only for qualitative visual inspection since they 

are less reliable as we discussed in Hu et al. [2014]. Several scalar quantities representing 
the total magnetic flux and magnetic helicity contents within certain boundary A = Ab 
are given in Table 1 together with the approximate average twist estimates thiTf- The 
other estimates for average twist ((r) and Tq) are based on more quantitatively reliable 
calculations to be described below. 

In our latest study of Hu et al. [2014], we performed systematic study of field-line twist 
distribution within ICME flux ropes based on the GS reconstruction method. The field- 
line twist, t{A), also as a single-variable function of A, is obtained by the graphic method 
described in Hu et al. [2014]. That is for each individual held line lying on a distinct 
surface of one particular A value, usually an open-ended cylindrical surface of closed side, 
denote the axial length along which the held line completes one full turn, in AU, then 
the field-line twist is simply 


r(A) 


1 


( 1 ) 
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in unit of turns/AU. The field-line length Lg is easily obtained by the line integral along 
each individual held line, i.e., by summing up the distances between the adjacent end 
points of each line segment. The results for the 7 events by this approach are shown 
in Figure 2. In the top panel, the mean value of each curve (each event) {t{A)) and 
the corresponding standard deviation At as given in Table 1 are shown as the square 
symbol and the associated error bar, respectively. They generally indicate that for most 
events, the twist remains fairly constant, excluding the region close to the center where the 
shifted hux function A' = |A — AqI = 0. Correspondingly, the length distribution shows a 
viable rate of increase with increasing A'. Note that each curve ends at certain value of 
A oi A' = Ac corresponding to the outermost closed loops on the flux-rope cross-section, 
represented by the equi-value contours of A. Beyond this boundary, the held lines can no 
longer complete one full turn within the computational box. Hence no held-line twist and 
length estimates by the graphic method are available [Hu et ai, 2014]. 

To circumvent this limitation and after observing that the twist distributions exhibit 
a trend of remaining fairly constant throughout the outer region of a hux rope, as hrst 
reported in Hu et al. [2014] and further demonstrated here, we employ a theoretical, 
constant-twist hux-rope model to assist in the analysis. To reinforce and justify this addi¬ 
tional approach, we put the results for all the events we have examined in Hu et al. [2014] 
and the present study together onto Figure 3, showing the average twist and associated 
standard deviations as they vary with Ac- The mean and median values of all points 
are 4.0 and 3.6 turns/AU, respectively. If the point of the largest standard deviation is 
excluded, they become 3.8 and 3.3 turns/AU, respectively. For half of the events of aver¬ 
age twist less than the median value, the standard deviations are small, indicating a hat 
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profile of t{A). Another general trend is that the larger size the flux rope is as indicated 
by larger value of A^, the smaller and less variable the twist becomes. 


2.2. Constant-Twist Gold-Hoyle (GH) Flux Rope Model 

The constant-twist or so-called Gold-Hoyle (GH) flux-rope model was originally de¬ 
veloped by Gold and Hoyle [I960]. It possesses rather simple and elegant forms for the 
magnetic held components in axi-symmetric cylindrical coordinate (r, 0, z)[Farrugia et ai, 


1999] 


B. 

B^ 


Bp 

1 + T2r2 
BpTpr 
1 + ■ 
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( 3 ) 
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Here the held-line twist by dehnition, = Tp = 2nTp, is strictly constant and is in 
the unit of radians/AU, which is also a signed quantity indicating the chirality of the 
hux rope. The parameter Bp corresponds to the axial magnetic held at the center of the 
hux rope (r = 0) which is set to be Bzp from the GS results as given in Table 1. They 
usually correspond to the maximum axial held during the interval (see Figure 1). The 
center of the hux rope is determined from the GS result as well and since we are only 
interested in deriving an approximation of held-line length as function of A, we don’t 
need to explicitly calculate r. The length can be expressed explicitly as a function of 
A thus can be directly estimated for each A value obtained from the GS reconstruction. 
The other parameter, tq, also given in Table 1 for each event, is obtained by taking the 
mean value of r for the outer loops where the twist variation is minimal, excluding the 
central core of each hux rope as we discussed earlier in Section 2.1. Largely based on the 
GS reconstruction results, the GH model provides an alternative and additional means 
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of estimating, especially extrapolating field-line lengths some of which are not available 
through the direct GS model estimates. 

From the GH model, because of the simple forms of the magnetic held components and 
the axisymmetric geometry, a hux function can be derived analytically 

A{r) = ~ 1.1(1 + ry). (4) 

Subsequently, the held line length per AU (i.e., for a section of the cylinder with an axial 
length Lgff = 1 AU) can be written as a function of A (Tq = 2nTo) 

Lgh = = ^l + ToV, (5) 

which tends to increase linearly with radial distance r from the center of the hux rope 
when Tor 3> 1. It is also worth noting that the GH model corresponds to a nonlinear 
force-free conhguration with the non-constant force-free parameter a = , varying 

with radial distance, i.e., along A shells as well, as originally derived by Gold and Hoyle 
[I960]. 

Table 2 summarizes the analysis results of measured and derived magnetic held-line 
lengths inside the selected MGs examined by Kahler et al. [2011a]. The entries of Date 
(1st column). Type III radio emission times (2nd column), measured held-line path lengths 
Lf. and D (3rd and 4th columns) are taken from Table 1 of Kahler et al. [2011a]. The 
path lengths D were obtained via the inverse-beta approach and were deemed inferior 
to the measurements Lg by the direct travel-time dispersion analysis. There are a few 
unacceptable values of D< 1 AU. The uncertainty in Le presented here is owing to the 
uncertainty in the exact timing of the arrival of the energetic electrons. The last two 
columns list the corresponding estimates of held-line lengths based on the direct GS 
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reconstruction output, Lg, and the GH model approximation, Lgh, respectively. The 
latter is obtained by applying the equation (5) with the necessary parameters supplied 
by the GS reconstruction results, i.e., the parameters Tq and Bq = from Table 1 and 
A = A'. The corresponding electron burst onset times at 1 AU are also given in the 
2nd column inside the parentheses. Note that only the event dates and times within the 
GS reconstruction intervals as indicated in Table 1 are listed for which our analysis can 
yield at least one estimate of Lg and Lgh- The others are considered outside of MG 
interval hence no analysis results are available. For the times listed which correspond 
to locations inside the MG but outside the loop boundary A' = A^ the estimates of Lg 
are not available while the estimates based on the GH model approximation can still be 
obtained. We defer detailed comparisons among these length estimates and discussion of 
their implications to Section 4. 

3. Case Studies 

In what follows three events are chosen to be presented as detailed case studies. The 
event 1 and 2 are selected because they possess the maximum number of electron burst 
onsets inside the MGs among all the events. The event 7 also contains a modest number 
of electron onset times and represents an extreme case of relatively and persistently long 
measured path lengths Lg throughout the MG interval. Thus these events facilitate a 
direct and broad comparison between measured Lg and estimated path lengths based on 
the GS reconstruction results and the GH model approximations. 

3.1. Event 1: 18 October 1995 


DRAFT 


February 19, 2015, 1:21am 


DRAFT 



X - 14 


HU ET AL.: MAGNETIC FIELD-LINE LENGTHS IN FLUX ROPES 


This event was also presented in Larson et al. [1997] and Kahler et al. [2011a], which 
possesses the maximum number of electron burst occurrences throughout the MC interval. 
The in-situ signatures of an MC structure are also strong, as seen from Figure 4a. The 
magnetic held magnitude is elevated and remains around 20 nT, the rotation in the GSE- 
Z component is the largest and clearly seen, and the plasma (3 is fairly low ~ 0.1, even 
after taking into account the electron temperature contribution {T^/Tp ~ 5). This is a 
relatively strong and long-duration MC event with a constant speed prohle and dominant 
magnetic held, indicating a typical hux-rope type magnetic structure embedded. This 
event was also examined by Hu and Sonnerup [2002] as a typical MC event to showcase 
the hrst application of the GS reconstruction method to the large-scale quasi-static MC 
hux-rope structures observed in-situ at 1 AU. The general presentation of the GS results 
is shown in Figure 4b and c. The data and a functional htting to the transverse pressure 
Pt = p + B'^/2fj,Q, the sum of the plasma pressure and the axial magnetic pressure, versus 
the hux function A as obtained along the spacecraft path at y = 0 are given in Figure 4b, 
together with a htting residual Rf indicating the goodness-of-ht of Pt{A) [Hu et al, 
2004]. Figure 4c shows the typical presentation of the cross-sectional map of the hux- 
rope structure as a contour plot of A{x, y) with the axial component superposed in 
color. It can also be viewed as a projection of the winding magnetic held lines lying on 
diherent iso-surfaces of A {A shells) onto the {x,y) plane. Figure 5 shows a rending of 
the 3D view with a few selected held-lines including the ones rooted on the locations of 
electron burst onset observations inside the surface A' = A^. Therefore the projected held 
lines that complete multiple turns around the axis will appear as the closed loops on 
the cross-sectional map of Figure 4c enclosed by the outermost loop highlighted in white 
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where A' = Ac. There are hve incidences of electron burst onsets along y = 0 as marked 
by cross signs with three occurring inside (in black) and the other two outside (in white) 
of the white loop. 

Figure 6 shows the double-axis plot of the distributions of held-line twist estimates (left 
axis) and the corresponding lengths (right axis) versus the shifted flux function A — Aq, 
including the available measurements of Lg with uncertainties, scattered at different A 
shells within the flux rope. Note that the shifted flux function is signed in this plot with 
the flux-rope center always located at A = Aq. Therefore the sign of the shifted flux 
function simply indicates the chirality of the flux rope: negative (positive) means right 
(left)-handed. The black thick curve and the three colored thin curves (see the legend in 
the top right-hand corner) represent the held-line twist estimates based on the graphic 
method and the other three approximate methods utilizing the magnetic hux and relative 
magnetic helicity content estimates, as described in details by Hu et al. [2014]. The 
graphic method yields the most accurate estimate but is limited to the inner region of 
loops satisfying A! < Ac. The results for the other three methods are only for reference 
purposes to visually inspect whether they follow the graphic method and the general 
trend of the twist distribution beyond the boundary where the graphic method ceases to 
provide an estimate [Hu et al, 2014]. As discussed earlier in Hu et al. [2014], the estimate 
by —d^p/d^t (green curve) would exhibit erroneous behavior of rapid rise toward the 
boundary of the flux rope (large A' values), as seen here, due to the rapid decrease in the 
estimate of d^ti but not in d^p. Overall, the twist distribution remains fairly low and 
constant with larger variations near the flux-rope center, yielding tq = 1.6 turns/AU for 
this case as indicated by the horizontal dashed line. 
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The corresponding length estimates (magenta curves and dotted curves) and measure¬ 
ment of Le (thick black horizontal and vertical lines) are overplotted versus the shifted 
flux function with the scales given by the right axis. There are two sets of estimates in 
this case. The thinner ones rise from 1 AU at A = Aq and increase toward the outer 
loops and they correspond to the estimates by using the default value L^ff = 1 AU. They 
do not intersect the measured Le at the locations along the A shells marked by the cross 
signs. For this particular case, there were five incidences of electron burst onset occurring 
within the GS interval with one occurring very close to the flux-rope center. There is 
a short vertical thick black line and a cross at A — Aq ~ 0, indicating the range and 
average of Le measured at that location. This measurement at this location enables us to 
determine Le// for this flux-rope event since the held line at the location near the center 
of the hux rope is mostly straight. Therefore a direct measurement of Le// is obtained 
in this case by taking the mid-point of the range of Le measured, Lg// ~ 1.3 AU. Then 
the actual held-line length estimates of both L* and Lgh are raised from their default 
values by simply multiplying the Lg// determined, resulting in the set of thick magenta 
and dotted curves. These corrected values will be used in the summary comparison with 
measurements. They now intercept a majority of vertical thick black lines except the one 
of the largest Lg at a location near the boundary. This adjustment by an ^eff > 1 AU 
only applies when such a direct measurement of Lg is available near the center of a flux 
rope. We choose the criterion for such locations along the A shells satisfying |A —Ao| < 10 
T-m. 

3.2. Event 2: 18 September 1997 
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This MC event has a very long duration, about two and a half days as seen in Figure 7a. 
The speed is fairly low, around 300-350 km/s during the GS interval. There are signihcant 
variations in the proton temperature Tp (the black trace in the third panel of Figure 7a) 
which does not show clear decrease inside the GS interval compared with the Tp values 
immediately outside, indicating the possible presence of signihcant plasma pressure. This 
results in a fairly modest plasma /9 ~0.5 within the GS interval, which is still depressed 
due to the relatively strong magnetic held magnitude. The plasma pressure becomes 
comparable to the axial magnetic pressure near the middle of the interval as shown in 
the last panel of Figure 7a, albeit there is less variation (smaller gradient) in the plasma 
pressure. The corresponding Pt{A) plot and the corresponding cross-sectional map are 
given in Figure 7b and c, respectively. Again this is a right-handed hux rope with a cross- 
sectional size of about 0.25 AU across. The inner loops enclosed by the white thick loop in 
Figure 7c occupies an area of a diameter about 0.1 AU. In this case, all the electron onset 
locations are outside of the closed loops bounded by the white loop where A' = Ac except 
for one point barely touching this boundary. Therefore most of the length estimates have 
to be obtained by the GH model-based extrapolation. 

The corresponding results including the held-line twist distribution and the actual mea¬ 
surements Lf. are shown in Figure 8, in the same format as Figure 6. The twist distribution 
remains fairly constant, especially in regions farther away from the flux-rope center, yield¬ 
ing To ~ 3.6 turns/AU. The held-line length estimates Lg rises from 1 AU at A' = 0 and 
increases to about 1.9 AU at A' = Ac, matching the measured Lg at that location. Beyond 
that point, no estimates of Lg are available, but the estimates by Lch are able to continue 
as illustrated by the dashed curve as A' increases toward the outer boundary of the hux 
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rope. These estimates seem to match the additional measurements of L^. except for the 
last point (left-most vertical bar) which is signihcantly lower than the estimated value 
Lgh ~ 3.8 AU, denoted by the cross sign at top. In this case, since there is no electron 
burst onset measurements close to the flux-rope center, the axial length of the flux rope is 
unknown and the default value Le// = 1 AU is used to obtain the corresponding held-line 
length estimates from the flux-rope models. The agreement with the measurements Le 
is reasonable. For most events examined in this study, we have to adopt this approach. 
Event 1 presented earlier and event 7, to be presented in the following subsection, are the 
only two exceptions. 

3.3. Event 7: 30 August 2004 

Event 7 is also a relatively large-scale event with a duration a little less than 24 hours, 
resulting in a relatively large-scale MC flux-rope structure. The in-situ data given in 
Figure 9a indicate a typical MC event: clear enhancement of the magnetic held magni¬ 
tude and rotation in direction, low proton temperature and low proton (3 within the GS 
interval. Although the magnetic pressure still dominates, because the ratio Te/Tp reaches 
10 in the GS interval, the plasma (3 is modest and in the range 0.1-1.0, owing largely to 
the contribution by the electron temperature to the total plasma pressure. The GS re¬ 
construction results including the contributions of both Te and Tp are shown in Figure 9b 
and c, in the same format as before. The Pt{A) curve shows a slight bend-over near the 
end to the right, which corresponds to the center of the flux rope as represented by the 
maximum A value. This behavior indicates a slight decrease in axial current density thus 
a weaker transverse held in the center. The corresponding cross-sectional map in panel (c) 
rehects this behavior with the transverse held nearly vanishing near the center whereas 
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the axial field Bz maintains a strong and fiat distribution inside a large area enclosed by 
the inner white loop. Such a configuration indicates that the field lines near the center are 
nearly straight with low twist. The 3D view of field lines is shown in Figure 10 where the 
inner field line, for example, the one of magenta color, is winding along the axis to large 
distance, ~ 1 AU, before completing one full turn. The two black lines rooted on two 
electron onset locations (two crosses inside the inner white loop on Figure 9c) near the 
center show similar behavior to the magenta line. Both the field-line length estimates Lg 
and Lqh for these two locations are available while only the estimate Lgh foL the other 
location outside of the loop A' = is available. Figure 11 shows the twist distribution 
and various length estimates along the A shells, similar to Figure 6. Here the adjusted 
length estimates are also given for an Te// ~ 3.0 AU, based on one measurement of Le 
close to the flux-rope center where A' < 10 T-m. The average twist is fairly low, yielding 
To ~ 1.2 turns/AU for the outer loops, the lowest among all the events. The one elec¬ 
tron burst onset measurement near the center yields a unusually large axial length of the 
flux rope, but the accordingly adjusted length estimates (thick magenta and dotted lines) 
show better agreement with measurements, matching 2 out of 3 values of L^. The one 
mis-match at the far left is almost outside of the flux-rope boundary defined hy A = Ai, 
beyond which the flux-rope interpretation based on the GS solution is less reliable. In 
other words, the location of this point could be outside of the MC flux rope and shall be 
excluded from the field-line length comparison. 

4. Summary and Interpretation of Results 

In this section, we summarize our analysis results presented in Table 2 and make direct 
comparison between the measured path lengths L^. (and D) and the derived ones from the 
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direct GS model output Ls and the constant-twist GH model estimates Lgh- For these 
handful of events, the path lengths obtained from the energetic electron burst onset mea¬ 
surements are in the range of 1 to 4 AU. A few exceptions exist for results corresponding 
to D which are less than 1 AU, thus deemed unacceptable. The apparent limitations and 
pitfalls of obtaining D based on energetic electron beams dispersion were discussed in sev¬ 
eral works [e.g., Kahler and Ragot, 2006; Wang et ai, 2011] but will not be repeated here. 
We adopt the results published by Kahler et al. [2011a] and their approach of weighing 
more the measurements of as better approximations of held-line path lengths. 

The derived path lengths from the GS together with the GH flux-rope models are within 
the same range as Le but are subject to an uncertainty in the effective length, Te//, the 
length of a section of the inhnite long cylinder that would correspond well to the flux-rope 
structure and the intrinsic characteristic quantities. Therefore the actual held-line length 
estimates are obtained by multiplying the lengths given for a section of unit axial length 
(usually 1 AU) by Le// in AU, whenever such a determination is available, as described 
in the case studies of events 1 and 7 presented in Sections 3.1 to 3.3. The uncertainty 
estimates in Lg and Lqh are based on errors propagated from the uncertainties associated 
with the measured electron onset times within the GS intervals. 

Figure 12 shows the ensemble distribution of measured held-line path length L^. along 
the A shells and the one-to-one comparison between Lg (and Lgh) and Le for all events. 
Figure 12a shows collectively all the measured Le along the A shells within GS intervals 
and their associated uncertainties. They exhibit a general trend of increasing path lengths 
with increasing A' = |A —AqI, i.e., with increasing radial distance away from the hux-rope 
center where A' = 0. For the events located near the center (to the left of the vertical 
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dashed line of A! = 10 T-m) the path length Lg would represent a direct measurement 
of -he//- For example, the one closest to A' = 0 at Lg ~ 1.3 AU corresponds to event 1 
presented in Section 3.1. The ones clustered around 3.0 AU correspond to event 7 
discussed in Section 3.3. Their locations on the cross-sectional maps of GS reconstruction 
results are close to the center of the flux ropes where the held-line twist values are small. 
A few guide lines are also drawn to further elucidate the trend and the coverage of the 
constant-twist GH model estimates. From the GS reconstruction results, an ensemble 
of held-line twists is obtained and presented in Table 1 and Figure 3, for example, from 
which a mean value of twist, as well as the minimum and maximum values are obtained. 
They are utilized to provide an estimate of coverage by the area bounded by the curves 
based on equation (5) varying along A shells for a given constant twist. In Figure 12a, 
the set of blue (red) curves corresponds to the length distribution along A shells based 
on the GH model for a constant twist of the minimum, mean, and maximum value from 
all held-line twist estimates, respectively, for L^ff = 1 AU (2 AU). In particular, the 
length variations for the mean twist values are drawn by dashed lines. Therefore it can 
be seen that the majority of the measurements falls within the region with the lower and 
upper bound provided by the GH model of the minimum and the maximum twist and for 
Leff G [1,2] AU. One exception is the measurements from event 7 as we discussed earlier 
which might be an indication that the ehective length could reach 3 AU in extreme cases. 

Figure 12b shows the direct comparison of versus Lg with associated uncertainties. 
Due to the limitation of the direct held-line length estimate from the GS reconstruction 
results, only 9 pairs of data points are available (the 5th column in Table 2). It shows 
good one-to-one correlation, especially considering that the correlation may be further 
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improved because the low points of low Lg values beneath the dashed diagonal line could 
be raised by a possible correction of being multiplied by an L^ff > 1 AU. The same 
comparison with expanded length estimates including additional GH model estimates is 
shown in Figure 12c where the additional pairs of Lqh and Le are marked by a cross and 
in black (the last column of Table 2). The correlation deteriorates compared with panel 
(b). However the number of outliers from the one-to-one line is few, especially counting 
only the ones above the dashed line, about 3, out of a total of 18, for the reason discussed 
above regarding panel (b). For completeness, we also show the same set of results and 
comparisons with replaced by D (4th column in Table 2) in Figure 13. The alternative 
length estimates D were provided by Kahler et al. [2011a] without error estimates and 
were not used in their comparison. The agreement of various model length estimates with 
the measurements seem to degrade compared with the previous hgure. For instance, the 
number of points above the dashed line increases in both panels (b) and (c). 

5. Conclusions and Discussion 

In conclusion, we have examined the flux-rope structures embedded within 7 MC events, 
in particular the held-line length and twist distributions, based on the GS reconstruction 
method and the constant-twist GH flux-rope model. We carry out direct comparison of 
held-line length estimates with the unique measurements of held-line path lengths obtained 
from timing observations of energetic electrons traveling along individual held lines from 
Sun to Earth. We limit our analysis to the same set of MG events reported by Kahler et al. 
[2011a] and employ their published measurements of to facilitate a highly comparative 
study but with diherent hux-rope models. Our conclusion, somewhat in contrary to 
Kahler’s, is that the hux-rope interpretation of the magnetic structures embedded within 
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MCs is largely consistent with the analysis of direct comparison between the modeled 
held-line length estimates and the direct measurements L^.. The correlations between L^. 
and Ls (and Lqh) are well established as seen in Figure 12 and the held-line length does 
exhibit a general trend of increasing from the hux-rope center. Such a trend as displayed 
in Figure 12a as a function of the shifted hux function has general implication for a hux- 
rope structure independent of specihc models. On the other hand, we agree with Kahler 
et al. [2011a] and others on that such a comparison provides additional evidence for the 
inconsistency of Lundquist model in characterizing the hux-rope structures observed in- 
situ at 1 AU. As we indicated in Hu et al. [2014], the magnetic held-line twist distribution 
within MC hux ropes often exhibits inconsistency with the Lundquist model, but better 
supports the GH model of a constant twist (see Figure 3 and associated descriptions). The 
present study further supports the hndings of such inconsistency and provides additional 
support for the GH model by the direct comparison of held-line length estimates with the 
corresponding measurements. It is also important to show that some electron burst onset 
observations are able to provide a direct measurement of the axial length of the section 
of a cylindrical hux rope, a critical parameter for the existing hux-rope models. Based 
on our analysis of a limited number of events, we argue that under most circumstances, 
such a constraint on the ehective axial length of a cylindrical hux rope is L^ff G [1, 2] AU, 
which has signihcant applications for the relevant studies of deriving and relating various 
physical quantities to their solar sources. 

It might not be hard to perceive why the comparison with the Lundquist model failed. 
Based on the Lundquist model, the held-line length would increase to inhnity at the 
boundary at which the axial held vanishes by dehnition. Therefore the Lundquist model 
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would yield large path lengths toward the outer loops of a flux rope. On the other 
hand, the GH model length would possess a more modest rate of increase from the center 
to the outer boundary of the flux rope, approximately linearly with r as indicated by 
equation (5), remaining hnite. Therefore the correlation between Le and Lgh is more 
favorable. As we discussed in Hu et al. [2014], the underlying theoretical consideration 
for advocating the GH model is that it describes space plasmas in non-linear force-free 
state which is well preserved from its origination from the Sun, propagation through 
the interplanetary space to reaching Earth. The ideal magnetohydrodynamic (MHD) 
conditions are probably well satished during the processes in the space plasma on the 
Sun and in the interplanetary space. The flux surfaces embedded within these structures 
remain distinct and well preserved upon their generation and are not destroyed by hnite 
and highly localized resistivity, resulting in a non-linear force-free state as observed in-situ 
at 1 AU. 

As an ongoing effort, we are extending the analysis to more events and utilizing more 
comprehensive sets of available observations. Some issues not addressed in the present 
work will be pursued in the forthcoming studies. For example, generally we would expect 
difficulty when the measured path lengths are exceedingly long and near the hux-rope 
center as we explained in the case study of event 7. Our interpretation of a hux-rope 
structure with an unusually long axial length of ~3 AU needs to be further validated by 
additional event studies. Another related issue is what ehect there is regarding the hnite 
plasma pressure gradient. A slight change in the model output of the conhguration of the 
hux rope would ahect the spatial locations of the electron onsets where the measurements 
of Le were taken. Such a change in location would yield change in the length estimates 
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by specific models. The amount of change may depend on whether or not the plasma 
pressure gradient is taken into account. Therefore a detailed assessment of the differences 
as resulted from the GS reconstruction results (non-force free in general) and the GH 
model (nonlinear force-free) estimates is planned for future work. 
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Figure 1. Summary plot of various physical quantities (unsigned) vs. the shifted flux function 
for the Wind MC events: (counterclock-wise from the top left panel) the poloidal (red pluses) 
and toroidal magnetic flux <hpi, the relative magnetic helicity, the held-line twist estimates th 
(red dots) and tf (blue dots) [Hu et al., 2014], the axial current, the axial current density, and 
the axial magnetic held. 


Table 1. GS Reconstruction Results of Selected Wind MC Events from Kahler et al. [2011a] 


Event # : GS Interval 
MM/DD/YYYY hh:mm:ss 

(r) ± At 
turns/AU 

Bzo 

nT 

lO^i Mx 

^p,max 

lO^i Mx 

K 

.''■r,max 

10^2 Mx^ 

Th 

Tf 

To 

1: 10/18/1995 18:59:30-(+l) 19:16:30 

1.57 ±0.26 

20 

1.2 

2.8 

2.9 

1.9 

2.2 

1.6 

2: 9/18/1997 03:55:30-(+2) 16:45:30 

4.55 ± 1.71 

12 

0.36 

2.1 

0.64 

5.0 

5.9 

3.6 

3: 11/6/2000 23:08:30-(+l) 18:46:30 

2.16 ±0.31 

25 

0.90 

2.0 

1.9 

2.3 

2.2 

2.2 

4: 7/10/2001 18:31:30-(+2) 07:47:30 

6.66 ±2.13 

7.8 

0.039 

0.35 

0.012 

8.3 

9.0 

6.0 

5: 10/1/2002 00:49:30-19:31:30 

4.20 ± 0.54 

26 

0.22 

0.99 

0.23 

4.7 

4.4 

5.0 

6: 7/24/2004 ll:56:30-(+l) 07:01:30 

2.11 ±0.83 

22 

1.6 

2.9 

3.5 

1.4 

1.8 

1.6 

7: 8/30/2004 01:4:30-21:42:30 

1.24 ±0.41 

12 

0.38 

0.55 

0.24 

1.7 

1.4 

1.2 
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Table 2. Measured and Derived Magnetic Field-line Lengths inside Magnetic Flnx Ropes'^ 


Date 

Type III (e-), UT 

Le, AU 

D, AU 

Ls, AU 

Lgh, AU 

18-Oct-1995 

19:56 (22:20:35) 

3.3-4.0 

3.30 


2.08 ±0.60 

19-Oct-1995 

05:18 (05:53:17) 

1.5-1.7 

1.06 

1.51-1.51 

1.50 ±0.003 

19-Oct-1995 

08:46 (09:30:58) 

1.2-1.6 

1.33 

1.40-1.47 

1.40 ±0.002 

19-Oct-1995 

10:28 (11:15:20) 

1.1-2.3 

1.65 

1.48-1.51 

1.48 ±0.02 

19-Oct-1995 

16:57 (17:25:50) 

1.7-2.2 

1.45 


2.16 ±0.06 

18-Sep-1997 

16:06 (16:52:37) 

2.6-3.2 

1.69 


2.63 ±0.05 

18-Sep-1997 

17:09 (17:56:00) 

2.7-3.1 

2.05 


2.38 ±0.04 

18-Sep-1997 

19:51 (20:19:57) 

1.8-2.1 

1.33 

1.89-1.94 

1.90 ±0.02 

20-Sep-1997 

03:16 (03:53:15) 

2.1-2.8 

1.33 


3.85 ±0.12 

07-NOV-2000 

00:08 (00:56:14) 

1.1-1.6 

1.08 


2.13 ±0.54 

07-NOV-2000 

15:40 (16:33:22) 

1.2-2.1 

0.98 

1.45-1.55 

1.41 ±0.04 

lO-Jul-2001 

22:53 (+1 00:44:47) 

1.5-2.5 

1.41 


1.85 ±0.26 

12-Jul-2001 

01:11 (02:28:27) 

2.7-3.6 

2.05 


1.53 ±0.02 

Ol-Oct-2002 

09:12 (11:55:21) 

1.7-2.7 

1.57 

1.11-1.35 

1.22 ±0.04 

24-Jul-2004 

18:43 (19:16:51) 

1.3-1.5 

1.14 

1.07-1.09 

1.15 ±0.01 

30-Aug-2004 

03:09 (03:44:47) 

1.9-3.0 

0.54 


3.49 ±0.10 

30-Aug-2004 

16:13 (17:47:05) 

3.0-3.4 

3.01 

3.13-3.26 

3.10 ±0.13 

30-Aug-2004 

18:09 (18:57:40) 

2.7-3.4 

3.31 

3.29-3.41 

3.17 ±0.12 


The Date, Type III times, and D are taken from Table 1, Kahler et al. [2011a]. 
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Figure 2. Summary plot of field-line twist r (top) and length (bottom) distributions (both for 
Leff = 1 AU) vs. the shifted flux function for the Wind MC events. Different colors represent 
different events as indicated by the legend in the top panel. The square symbol and associated 
error bar at the end of each curve where A' = A^. indicate the mean and the standard deviation 
of each curve. 
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Figure 3. The mean and standard deviation of held-line twist t{A) of 25 magnetic flux-rope 
events. Each data point with associated error bar is plotted at the corresponding value of Ac. 
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FigurG 4. The GS reconstruction result for event 1 in Table 1. (a) Time series of Wind spacecraft measurements: 
(from top to bottom panels) the in-situ magnetic field magnitude (black) and GSE-X (red), Y (green), and Z (blue) 
components, the plasma bulk flow speed, the proton density (left axis; blue) and proton (black) and electron (green; if 
available) temperature (right axis), the plasma /3 (black) and the electron over proton temperature ratio (red; if available), 
and the plasma and axial magnetic field (red) pressure. The vertical lines mark the GS reconstruction interval as given 
beneath the last panel, (b) The measurements of Pt{x, 0) versus A{x, 0) and the fitted Pt[A) curve (thick black line). The 
flux rope boundary is marked a,t A = Ab and a fitting residue Rf is denoted, (c) The cross-sectional map of the solution 
A{x,y) (black contour lines) and the axial field Bz{A) (filled contours in color). The yellow arrows are the measured 
transverse magnetic field along the spacecraft path (t/ = 0). The white contour lines denote the boundary A = Ab (outer) 
and A! = Ac (inner) while the white dot denotes the center where the axial field is the maximum and A = Aq. The crosses 
along y — denote the locations where the electron burst onsets were observed. The ones inside the white loop are in black. 
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zs=(0.4642 -0.756 -0.4614) 



Figure 5. 3D view of the flux-rope structure toward Sun for Event 1 with selected held lines. 
The view angle is such that north is upward and the ecliptic plane is horizontal. Black lines are 
the held lines rooted at the footpoints where the electron onsets were observed. Circles mark the 
locations where the held lines complete one full turn around the z axis. The orientation of z axis 
is given on top in GSE coordinate. 
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x„=1.57 



Figure 6. Field-line twist and length distributions along the A shells for Event 1. The twist 
values are indicated by the left axis with the estimate by the graphic method in thick black curve. 
The results from the other three approximate methods are given by the thin blue, green and red 
curves, respectively, as indicated by the legend. Each curve ends at A = Aj,. The horizontal 
dashed line denotes tq as given on top of the plot in turns/AU. The magenta lines represent the 
length estimate Lg, the dotted lines Lgh, and the black thick verticle lines Le with uncertainties. 
All the length scales are given by the right axis. The crosses mark the locations along the A 
shells where the measurements of Lf, were obtained, and the corresponding estimates of Lgh- 
The set of thinner magenta and dotted curves originating from 1 at A' = 0 corresponds to the 
default value L^ff = 1 AU. The other set corresponds to Lg// = 1.3 AU. 
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Figure 7. The GS reconstruction result for event 2 in Table 1. Format is the same as Figure 4. 


X|j=3.58 



Figure 8. Field-line twist and length distributions along the A shells for Event 2. Format is 
the same as Figure 6. Here the effective axial length takes the default value Tg// = 1 AU, only. 
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Figure 9. The GS reconstruction result for event 7 in Table 1. Format is the same as Figure 4. 
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zs=(0.4343 0.5224 -0.7338) 



Figure 10. 3D view of the flux-rope structure toward Sun for Event 7. Black lines are the 
held lines rooted at the footpoints where the electron onsets were observed. Format is the same 
as Figure 5. 



Figure 11. Field-line twist and length distributions along the A shells for Event 7. Format 
is the same as Figure 6. Here the set of thicker magenta and dotted curves corresponds to 
Leff ~ 3.0 AU. 
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(a) 



(b) (c) 



Figure 12. Summary and comparison of field-line length estimates with Le'. (a) the ensemble 
of measurements vs. A' for all events; the dashed and dotted lines mark the variations of 
held-line lengths of GH model with certain constant twist, and the vertical dashed line denotes 
A' = 10 T-m (see text for details), (b) the one-to-one plot of Lg (in magenta) vs. Le, and (c) the 
one-to-one plot of both Lg and Lgh vs. Le where the latter is marked by black cross signs and 
vertical lines. In panels (b) and (c), the dashed line denotes the one-to-one diagonal line while 
the dotted lines mark a 10% uncertainty bound. 
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(a) 



(b) (c) 



Figure 13. Summary and comparison of field-line length estimates with D which were given 
by Kahler et al. [2011a] without associated uncertainties. Here a uniform 10% uncertainty in D 
is assumed. Format is the same as Figure 12. 
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